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PRECISION GPS EPHEMERIDES AND BASELINES 
NASA GRANT NO. NAGW-2717 


Final Report 

January 1, 1991 - June 30, 1992 


INTRODUCTION 

The emphasis of this grant was focused on precision ephemerides for the Global Positioning 
System (GPS) satellites for geodynamics applications. During the period of this grant, major 
activities were in the areas of thermal force modeling, numerical integration accuracy improvement 
for eclipsing satellites, analysis of GIG ’91 campaign data and the Southwest Pacific campaign data 
analysis. 

THERMAL FORCE AND ECLIPSING ORBITS 

Papers on thermal imbalance force modeling and eclipsing orbit analysis were presented at the 
AIAA Space Flight Mechanics Conference held in February 1992 at Houston, Texas. Details are 
included in Appendix VII and Appendix VIII of the Final Report for NASA Grant No. NAG5-940. 

GIG ’91 DATA ANALYSIS 

An intensive global GPS observation campaign called "The First GPS IERS and Geodynamics 
Experiment - GIG ’91" was organized by the Jet Propulsion Laboratory for a period of three weeks 
during early 1992. The campaign spanned part of GPS Weeks 576, 577, 578 and part of Week 579 
starting from January 22, 1992, to February 13, 1992. About 120 stations distributed globally, 
collected data during this period using a variety of GPS receivers which included Rogues, TI-4100’s, 
Trimbles, Ashtechs, Mini Mac 2816 AT’s, WM-102’s; in addition, data from the five continuously 
operating DMA stations were also available. From among these, data from 20 globally distributed 
stations were chosen for orbit determination and Earth orientation parameter determination 
experiments at UT/CSR. Of these, 17 sites had Rogue receivers, two TI-4100 and one Mini Mac 
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2816AT. Orbit determination experiments included computation of orbits using data from a reduced 
set of global stations and from an expanded set, and comparing the baselines and other criteria to 
assess the orbit accuracy. In addition to the GIG’91 data set mentioned above, data collected during 
the 1989 South West Pacific Campaign also were used. Results of this study were presented at the 
1991 AGU Spring Meeting held at Baltimore, Maryland, during May 1991. Further details are 
included in the Appendix. 

Using the same GIG’91 data set several experiments were performed at CSR to determine the 
pole positions (X p , Y p ) during this time period. The main objective of the experiments was to 
explore strategies for determining the Earth orientation parameters (EOP), polar motion and UT1, 
using global GPS data. Short arc and long arc with sub arc parameters were considered as possible 
approaches. The estimated pole positions were compared with SLR and VLBI determined polar 
motion series in order to assess the quality of these determinations. Results of this study were 
presented in the special workshop held in Ahrweiler, Germany, during August 1991 and at the XX 
General Assembly of the IUGG during August 1991 in Vienna, Austria. A copy of the presentation 
is included in the Appendix. 

SWP CAMPAIGN DATA ANALYSIS 

Analysis of data collected during the South West Pacific (SWP) campaigns contributed to 
accuracy assessment of the GPS satellite orbits. The SWP campaign data were collected during the 
summers of 1988, 1989 and 1990 and these data have been analyzed in various stages at CSR. These 
data sets complement the CIGNET data in forming a better global distribution and facilitates various 
orbit and baseline experiments. Studying the repeatability of baselines between various sites in the 
SWP network provides one method to evaluate the accuracy level of the computed orbits. In depth 
analysis of few weeks of data collected at the sites in Tonga region (SWP) during the above three 
years were performed at CSR and the results indicate that the computed GPS orbits yield baseline 
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repeatability at the 10 to 20 parts per billion (ppb) level. A brief summary of these results were 
presented at the 1991 AGU Fall Meeting held in San Francisco during December 1991. Copy of the 
presentation is included in Appendix. 

CONCLUSION 

Based on the research efforts at CSR in the area of precise ephemerides for GPS satellites, the 
following observations can be made pertaining to the status and future work needed regarding orbit 
accuracy. There are several aspects which need to be addressed in discussing determination of 
precise orbits, such as force models, kinematic models, measurement models, data 
reduction/estimation methods etc. Although each one of these aspects has been studied at CSR in 
research efforts under this (and the previous) grant, only points pertaining to the force modeling 
aspect are addressed here. 

Dynamic Modeling - Current Status 

At present the following known forces are modeled in routine computation of GPS satellite 
orbits: 

(1) Nonspherical Earth gravity acceleration represented by one of the extant gravity fields such 
as GEM-L2, GEM-T1 or GEM-T2 truncated to 8x8 is adequate. However, increasing the degree 
and order and/or tuning the field (at least the resonance coefficients) for GPS orbits may slightly 
improve the accuracy of this perturbation modeling. 

(2) represented as point masses are adequate. 

(3) force modeling and the geometric tide effects must be considered in the measurement 
models. 

(4) Perturbation due to solar radiation pressure is adequately modeled by the ROCK4 models; 
however, it is necessary to scale these accelerations by at least one adjustable parameter. 
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(5) Perturbation due to venting of the heat source in the GPS satellite is modeled as ’Y-bias 
acceleration’ scaled by an adjustable parameter. 

Values of weekly estimates of this scale parameter considered over a period of time do not 
exhibit any systematic trend except for a somewhat weak correlation with the eclipsing period. 
Experience shows that it is necessary to include this perturbation along with an adjustable parameter 
in order to obtain good fit of the data. However, in a long arc solution, considering several sub arcs 
for this perturbation does not significantly improve the rms of fit 

Following are some of the factors which could contribute to additional improvement of orbit 
accuracy beyond the current level. 

(1) Perturbation due to imbalance in thermal radiation has been shown to cause differences in 
orbit prediction at the level of a few meters over a period of about one week or more. Although 
small, inclusion of this perturbation may help in achieving baseline accuracies (and/or 
repeatabilities) at parts per billion level. But the difficulty in considering this perturbation routinely 
is due to the fact that nonlinear partial differential equations (heat equations) must be solved 
simultaneously with the ordinary differential equations of motion, which causes significant 
complications in algorithm and computation even for a modest approximation of the satellite 
configuration. Hence, careful evaluation of the costs and benefits of including this perturbation is 
needed. 

(2) handling discontinuities in function values (occurring in SRP acceleration at shadow 
crossings), can be overcome by the ad hoc modification of the integrator back difference table. 
However, inclusion of this modification did not seem to improve the rms of fit or the prediction error 
in real data processing, although improvements were obvious in simulation studies. The reason for 
this anomaly is not known at present and will have to be investigated before this feature can 
routinely be included in orbit computation. 
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Dynamic Modeling - Future Study 

There are indications (evidenced by discontinuities in daily/weekly solutions and by prediction 
errors) to the effect that all the perturbations described above do not completely or exactly represent 
all the forces acting on the GPS satellites. There may be other unmodeled forces such as 
unintentional thrusting (due to outgassing, momentum dumping, attitude correction etc.) or due to 
other natural phenomena. One of the ways in which such unknown and unmodeled perturbations 
could be accounted for in orbit computation is to estimate empirical accelerations. Such an approach 
needs detailed analysis in the future. 
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APPENDIX 


1. "GPS Ephemeris Accuracy Improvement From Global Data Set," AGU 1991 Spring 
Meeting, May 28-31, 1991, Baltimore, Maryland. 

2. "Analysis of GIG 91," presented at GIG Workshop, Ahrweiler, Germany, August 1991. 

3. "Global GPS Orbit Determination," XX General Assembly, IUGG, August 1 1-24, 1991, 
Vienna, Austria. 

4. "Geodetic Analysis of GPS Measurements Near the Tonga Trench: 1988-1990," AGU 
1991 Fall Meeting, December 1991, San Francisco, California. 

5. "GPS Reference Frames and Earth Rotation," Proc. Sixth International Geodetic Symp. 
on Satellite Positioning,, Columbus, Ohio, March 17-20, 1992. 

6. "GPS Orbit Accuracy," Proc. Sixth International Geodetic Symp. on Satellite 
Positioning,, Columbus, Ohio, March 17-20, 1992. 

7. "GPS Reference Frames and Orbit Accuracy," AGU Spring Meeting, Montreal, Canada, 
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GPS EPHEMERIS ACCURACY IMPROVEMENT 
FROM GLOBAL DATA SET 
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GLOBAL TRACKING NETWORK DATA FOR WEEK 499 
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SOLUTION CHARACTERISTICS FOR THE GLOBAL DATA SET 
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TROMSO-ONSALA VECTOR BASELINE ACCURACY 
COMPARISON WITH VLBI SOLUTION 
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RMS OF TRAJECTORY DIFFERENCES 
9 GS vs 20 GS SOLUTIONS 
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ECF EPH DIFFERENCES (WEEK 499) 
(UT/CSR DD PHASE SOLN - NSWC EPH) 
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GIG' 91 Workshop 

Agenda 


Monday , 5 August 1991 

Arrival of the participants 

19.00 - 20.00 Open-air reception by good weather otherwise in the restau- 
rant of the hotel 


Tuesday, 6 August 1991 

08.00 - 09.00 Registration 

09.00 - 10.30 Opening Session (Chair: H. Seeger, IFAG) 

Greeting and Remarks from Hosts: H. Seeger (IFAG) 

Greeting and Remarks from IERS-Background and objectives of 
the IERS GPS Campaign: M. Feissel? (IGN) 

Greeting and Remarks from the IERS GPS Technique 
Coordinator: W. Melbourne (JPL) 

IERS Latest UTPM Results for the GIG Experiment Period: 

M. Feissel? (IGN) 

10.30 - 11.00 Coffee break 

11.00 - 12.30 The Field Campaign (Ch: W. Schluter, IFAG/R. Neilan, JPL) 

Summary of the Campaign-Planning, Standards, Procedures, 
Operations, Performance, etc: R. Neilan/S. Fisher (JPL) 

Reports from Various Operations Teams (Representatives from 
organizations that participated in field ops will be asked to 
give a brief account of their experiences, insights, etc.) 

12.30 - 14.00 Lunch 

14.00 - 15.30 The Field Campaign (cont.) 

Continued Reports from Ops Teams 

Discussion: Lessons Learned, Recommendations for Future 
Campaigns (Standards, Procedures, Documentation, etc.) 

15.30 - 16.00 Coffee break 

16.00 - 18.00 Data Analysis - I: Preparation 

Status of Data Pre-Processing & Distrib., Site Ties 
(S. Fisher) 

Reports by Pre-Processing Centers (All) 

Summary of Analysis Standards Proposed by JPL (G. Blewitt) 
Discussion of above and related topics 

Open-air grillparty 


19.00 



/ Uti.JU - 1U.JU 

h 

uara Analysis - 11 : Kesuius - ine plan is to organize this all-day session by ^ 

topics. Those with results in several areas will therefore give several short 
presentations. The main topics are: 


1. Descriptions of software and modeling strategies 

2. Baseline Solutions 

3. GPS Orbit Solutions 

4. Earth Orientation Solutions 

5. Geocenter Solutions 

6. Discussion, Comparison & Analysis of Results 


As of 25 July, the known groups planning to present results are: 
U. Texas - B. Schutz 
DGFI - K. Kaniuth 
NSWC - E. Swift 
UNAVC0 - C. Rocken 
MIT - R. King 

JPL - Blewitt, Lichten, Lindqwister, Webb, Yunck 
There is still space for more 

10.30 - 11.00 

Coffee break 

11.00 - 12.30 

Data Analysis - II: Results (cent.) 
(See Above) 

12.30 - 14.00 

Lunch 

14.00 - 15.30 

Data Analysis - II: Results (cent.) 

(See Above) 

15.30 - 16.00 

Coffee break 

16.00 - 18.00 

Data Analysis - II: Results (cont.) 
(See Above) 

19.30 

Bus departure for wine-tasting 

20.00 - 22.30 

Wine-tasting in Ahrweiler 

22.30 

Busses return to hotel 

Thursday, 8 August 1991 

08.30 - 10.00. 

Special Topics (Preliminary list): 


Multipath studies with GIG data (C. Rocken) 

Do high latitude sites pose special problems? (TBD) 
Implications of SA/AS for global GPS measurements (T. Yunck) 
Status of 3CAR GPS geodesy project in Antarctica (J. Manning) 
Surma ry of early GIG results (TBD) 

10.00 - 10.30 

Coffee break 

10.30 - 12.30 

Open Discussion (Suggested topics): 

v • *■■■■ 

Inplications of early GIG results 

Establishing a GPS-based global reference frame Integrating GPS into IERS 
operations and products Plans for a follow-up 616*91 vorkshop before IGS 
campaign? 

12.30 : 

Lunch 
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ANALYSIS OF GIG 91 
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Adjust all stations with 1 meter a priori, all other estimated parameters 
have no a priori 

Pole position error estimate (RMS): 1 mas 



GIG 91 DATA 
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GLOBAL GPS ORBIT DETERMINATION 
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SOFTWARE VALIDATION 
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SOFTWARE VALIDATION (continued) 
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GEODETIC ANALYSIS OF GPS 



DECEMBER 1991 
SAN FRANCISCO 




HI continental crust A A A subduction rone _TL_n spreading ridge offset by transform faults 

‘•••volcanoes uncertain plate boundary collision lone 

earthquake zone * movement of plate 21 Plate boundaries and active zones of the Earth's crust 














SOUTHWEST PACIFIC DATA CHARACTERISTICS 
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SOUTHWEST PACIFIC PROCESSING STATUS 
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Figure 2. 1988 Southwest Pacific Campaign Global Sites 



SWP-90 ANALYSIS STRATEGY AND EXPERIMENTS 
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SWP-90 EVALUATION SOLUTIONS 
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SWP-88 SOLUTIONS 


3 On 


wo .S 

Zt co 
cd 

M V 

1) co 

o irJ 
P* <D 

T3 S 

c c3 

* a 

3 ^ 
3 .a 

^ o 
£ T) 

£ & 
S 'S 

0 § 

1 . +- 1 > 

^ *i—i 

<2-e 

c/5 O 

c3 o 

a « 

cd a 

Sfi 

P-l CO 
<D ^ 
co 

TT 
c/i «* 

O c/i 

S3 § 

S*§ 


CO *— ■ i 

<N c4 
II II 
00 00 


’ I — ’ N 

<n O' ' 
21 CO CO 
o o Tt 
cd 


CO cd 

T3 o 

* fci 

E o 


OOr^ S 

<d p n 

ISK 

00 • • 


G M 
O o 


epeatabilities (combined 444 anc 

• Rarotonga/Tongatapu: 4.6 cm 

• Rarotonga/Vava’u: 4.1 cm 

• Rarotonga/W. Samoa: 6.9 cm 



CONCLUSIONS 


-a *ss 

<1> .5-h 
o q 


.2 || 

C/3 ^ 

g »P o3 

j.S's 

(2 c3 £ 


y 2 « J 
q H g g 
§ g £ <D 
*-• Jft *P <D 


+- > »ri 

13 g S a 

t! d c 3 o 

:< § 

a g> 8 §> 

g C/5 8 g 

S ^ - s 

•5 03 w ^ 

<D o o 13 

IK'S 

3 &-S S 

W &.T! ^ 


P o 

p £ 

W q 


o 00 


> 00 
5 o\ 


j'O o 

3 § 13 3 

5 g-81 

, £ -e «> S 

c d 8 O 

O « o w 

-g < * -y § 

3 x: « o 2 


3 xi w 

g o 9 

O £ 8 

u c* 


co o 2 

P5 £ S 

O -m ‘rt 
00 <D CO 

q C d> 


T) <D 
<L> > 

£>J3 

£ £* 

Cj r^J 

c3 is 

<D td 
oo cd 
^ <U 

g (X 

> P 

w> q 

ft 

•P w 
oo <_> 
P 03 

CO CO 

1=5 <0 
C O 

§ 3 

£ g 
*a j3 
<D fa 
cu-tj 
X T3 

w 


£ X) 

0 a, 

m Ph 
«N O 

1 1-1 

* Irt 

co cd 

P « 
£ g 

° 8 

10 P 

.t & 

> 60 
w 03 
_ C/3 

s I 

2 is 

^ CO 


>, o 

- 4 -* 


rP JP 
cd o 

td N 

£S 

•S3 *«-i 

•a 13 

CO 

ft 03 

•S»3 


co ft 

4=5 O 
ft co 
co *rj 


00 ft 
OO ft 

as o 


oo • 


oo • 


(Baseline: 2296 km) 




Proc . 6th International Geodetic Syrnp. on Satellite Positioning, Vol /, 
Ohio State Univ., Columbus, OH, March 17-20, 1992 


GPS ORBIT ACCURACY 


P. A. M. Abusali, B. Schutz, D. Kuang and H. Rim 
Center for Space Research 
The University of Texas at Austin 
Austin, TX 78712-1085 (USA) 


ABSTRACT 

Previous analysis of GPS data collected by special campaigns (e.g., CASA UNO, GIG-91, 
etc.) has shown that the satellites exhibit somewhat different characteristics, especially between the 
eclipsing and non-eclipsing satellites. Using the GIG-91 data which provided a reasonable global 
distribution of stations, the influence of unmodeled orbital effects has been examined using double 
differenced earner phase data. While there is evidence that suggests the unmodeled orbital effects 
are not a limiting factor in achieving a part in 10** level in baseline results, these components may 
be factors in reaching a part in 10**. Experiments with the GIG data set include comparison with 
other ephemerides and the estimation of empirical parameters for die purpose of improving the 
model error characterization. 

1. INTRODUCTION 

Previous results have shown that unmodeled forces exist when the GPS satellites are in 
eclipse season, i.e., the period during which the satellite experiences the umbra/penumbra of either 
the Earth or the Moon [Schutz et al., 1991; Fliegel et al., 1992; Gouldman et al. f 1989]. Possible 
contributors to the observed effects include the proper representation of the discontinuity associat- 
ed with the shadow boundary and the implementation of appropriate adjustments in the numerical 
integration algorithm [Lundberg et aL, 1991]. However, studies of this effect by Feulner et 
al.[1990] demonstrate that, while the effect can be significant, it does not account for most of the 
observed effect. Another effect that is associated with thermal radiation imbalance was examined 
by Vigue et al.[199 1] who demonstrated that the effect should be observable. 

The objective of this investigation was to assess the GPS orbit accuracy and to examine pos- 
sible parameterization to account for observed mismodeling of the measurements. Such an exam- 
ination cannot be accomplished with a regional network, and is best suited for a global tracking 
network. The global data set of the GIG-91 Campaign [Melbourne, 1992] offers an opportunity to 
examine a variety of aspects concerning the fidelity of the GPS force, kinematic and measurement 
models. The campaign used about 20 Rogue P-code receivers plus several TI-4100 P-code receiv- 
ers and numerous codeless receivers. 

2. DATA AND MODELS 

The software used for the data analysis was the set of programs know collectively as TEX- 
GAP, described by Schutz et al.[1992]. For the results of this paper, ionospherically corrected 



phase measurements were used in a double difference mode. 


The specific receivers used in the analysis were 17 Rogue receivers plus TI-4100 receivers 
in the Pacific to provide improved global coverage, as given in Table 1 . The GPS force models fol- 
lowed the current IERS Standards [McCarthy, 1989] and scale factors on the ROCK4 (Fliegel, 
1992] and y-bias parameters were estimated. The Chao [1974] troposphere model was used. La- 
geos-derived polar motion and UT1 were used as a priori. Pseudo-range measurements were used 
to verify and/or correct the respective receiver clocks. The reference frame is given by Schutz et 
al.[1992]. 

For this study, data from GPS Weeks 578 and 579 (days 34 to day 41, inclusive) were used. 
All available satellites (5 Block I and 10 Block II) were included in the analysis. 

3. ESTIMATION STRATEGIES 

All results were obtained using multi-satellite versions of UTOPIA, known as MSODP. The 
estimation process is based on a batch algorithm, using Givens rotations to solve the least squares 
problem. For the results, three sites were fixed (Goldstone, Wettzell, and Hobart) and the remain- 
ing 17 sites were adjusted. For each arc, satellite position and velocity were estimated at the initial 
time point of the arc, a solar radiation pressure scale factor and a y-bias were estimated. Zenith de- 
lay parameters were estimated at 2.5 hour intervals and phase ambiguities were estimated on each 
pass. In all cases, the a priori covariance was assumed to be infinite, thus allowing all parameters 
to freely adjust. 

Alternate empirical forces were introduced for some cases. These forces include radial, 
along-track and cross-track components represented by a periodic function. The period of this 
function was adopted to be the orbital period, thus the empirical force accommodated once/revo- 
lution effects. The estimated parameters were amplitude and phase of the function. 

The arc lengths included a series of “short arcs" of one-day duration, each of which was in- 
dependent of the other arcs. For this study, a “long arc” consisted of a 5-day arc in which a single 
set of orbital parameters for each satellite were estimated. For the one-day arcs, three cases have 
been examined: 

• Case 1 : estimate ROCK4 scale parameter and y-bias for each satellite 

• Case 2: estimate coefficients of once/revolution radial, transverse and normal perturba- 
tions instead of radiation pressure and y-bias 

• Case 3: same as Case 2 except a constant along-track perturbation was estimated instead 
of the once/revolution transverse force 

For the 5-day arc, a strategy similar to Case 1 was followed except that two y-bias parameters for 
each satellite were estimated. 

4. RESULTS 

The statistics of the Case 1 results are shown in Table 2. In general, the RMS of the double 



difference residuals from the one -day arcs were in the range of 1.2 to 1.6 cm. Examination of the 
raw phase measurements suggests that the ionosphere connected phase measurement has a preci- 
sion of about 0.3 to 0.4 cm, thus leading to the conclusion that the precision of the double differ- 
ence (DD) measurement should be about 0.6 to 0.8 cm. The discrepancy between the DD precision 
estimate and the values in Table 2 is indicative of one aspect of mismodeling. However, it cannot 
be concluded that the discrepancy is caused completely by orbit mismodeling and the possibility 
that measurement systematics, such as multi-path, are a contributor must be considered. 

The increased DD RMS from the 5-day arc, however, is indicative of a level of orbit model 
error since the measurement systematics are not dependent on the arc length, but force models are 
significantly dependent on the arc. Nevertheless, although the 5-day arc should use daily sub-arc 
values of SRP and y-bias to more nearly match the one-day arcs, past experience has shown that 
such representations do not substantially reduce the RMS on the long arc [Schutz et ah, 1990]. 

There are possible sources of orbital mismodeling: gravitational and nongravitational. Expe- 
rience with satellite laser range (SLR) measurements to the Etaion satellites, however, suggests 
that no significant gravitational mismodeling exists [Eanes, ct al., 1991]. The two Etaion satellites 
were launched into GLONASS-like orbits by the USSR in 1989. Both are spherical, with a reason- 
ably low area to mass ratio. The dominant model error on the Etaion satellites is nongravitational 
in origin, however, the nature of the nongravitational effects on Etaion is quite different than GPS 
and the Etaion experience cannot be extrapolated to GPS (or GLONASS). Concerning the gravi- 
tational contributions, the fact that the GPS satellites are in “deep resonance” distinguishes them 
from the Etaion satellites which are not; thus, there is still the possibility of a gravitational effect, 
but it is most likely of very long period and would not be evident in arcs with a duration of several 
days. 

The mismodeling is further evidenced by discontinuities in the common time point between 
the one-day arcs. For PRN-3, the differences between the Case 1 one-day arcs and the 5-day arcs 
is shown in Fig. 1 for the radial, along-track and cross-track components. As shown, the disconti- 
nuities at the common time point are several meters in some cases, while others are at the level of 
2 meters. The discontinuities are associated with mismodeling on the one-day arcs, however, the 
magnitude of the discontinuity is, in part, determined by the mismodeling on the 5-day arc. Note 
that the magnitudes of discontinuities in the radial and cross-track directions are much smaller than 
those in the along-track direction. 

For PRN-3, the 5-day arc was compared with the ephemeris produced by Defense Mapping 
Agency. The differences are shown in Fig. 2. This comparison was accomplished without any ad- 
justment to either the DMA or the UT ephemerides, and was formed by directly differencing the 
two ephemerides in the Earth-fixed system and transforming the difference into radial, along-track 
and cross-track components. Because of the difference in GM used in the ephemerides (DMA: 
398600.5; UT: 398600.441 knrVs 2 ), a radial bias exists at the meter level. The periodic differences 
probably reflect model differences, including reference frame differences. Since the two cases 
were generated by independent software and different global tracking networks as well as different 
data types, the differences can be regarded as an indication of the level of GPS orbit accuracy. The 
RMS of differences are 1 .6 m radial, 2.3 m along-track and 2.6 m cross-track. The size of the cross- 
track in comparison to the along-track is an indication of reference frame differences. 



Additional experiments using the once/revolution force model characterizations were con- 
ducted. The Case 3 result for PRN-3 is shown in Fig. 3 for the coefficients of the radial and cross- 
track components and the constant (over one day) along-track component is shown in Fig. 4, in- 
cluding the formal error of the respective daily estimates. Although the trends exhibited by these 
parameters appear to be systematic, it should be noted that the effect on the RMS of the one-day 

arcs has been small. Further experiments will be conducted using these parameterizations in long 
arcs. 

5. CONCLUSIONS 

Based on analysis of the GIG-91 data set, double difference phase residual RMS at the 1 .2 to 
1.6 cm level have been obtained for one-day arcs, while a 5-day arc shows 3.8 cm. The one-day 
arcs are probably influenced by both unmodeled forces on the GPS satellites and by systematic 
measurement model errors, while the 5-day arc is expected to be dominated by force model errors. 
Experience with other satellites at similar altitudes suggests the dominant force model error has a 
nongravitational origin. Comparison of the one-day arcs with the five-day arc shows discontinui- 
ties at the common time point of the one-day arcs with differences of several meters. Direct com- 
parisons with DMA ephemerides show differences at the 2 to 3 meter level (RMS), thus providing 
an indication of the orbit accuracy over days 34-38. Use of empirical force models as a means of 
investigating the nature of possible model errors was applied to one-day arcs, with results that ex- 
hibit systematic characteristics. Future studies will investigate these parameterizations. 
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TABLE 1. GIG-91 SITES 


Rogue Receivers: 

Yaragadee, Australia 
Canberra, Australia 
Santiago, Chile 
Hartebeesthoek, S. Africa 
Kokee Park, Hawaii, 

Usuda, Japan 
Goldstone, CA 
Victoria, BC 
Fairbanks, Alaska 

Minimac 28 16 AT Receiver: 

Hobart, Tasmania, Australia 

TI-4100 Receivers: 

W. Samoa 
Easter Island 


Wettzell, Germany 
Madrid, Spain 
Matera, Italy 
Kootwijk, Netherlands 
Ny Alesund, Norway 
Tromso, Norway 
Algonquin, Ontario 
Yellowknife, NWT 


TABLE 2. ARC STATISTICS 


One-Day Arcs 


Day 

Passes 

DD Observations 

RMS (cm) 

34 

221 

11545 

L550 

35 

224 

11421 

L511 

36 

291 

15035 

1.516 

37 

287 

14823 

1.396 

38 

280 

13792 

1.377 

39 

289 

12397 

1.378 

40 

304 

13383 

1.265 

41 

312 

13742 

1.204 

Five-Day Arc 

34-38 

1300 

6 6629 

3.857 


DD denotes Double Difference 



DIFFERENCE (M) DIFFERENCE (M) 



TIME (DAY) 


Figure 2. Difference between UT/CSR 5-day arc and DMA ephemeris of PRN03; 
Epoch: Feb.3, 1991. 





Figure 3. Estimated coefficients of once/rev radial and cross-trade empirical 

acceleration for eight daily arcs of PRN03. Day No.l is Feb.3, 1991. 



Figure 4. Estimated constant along-track empirical accderation (Case 3) for eight 
daily arcs of PRN03. Day No.l is Feb.3, 1991. 
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ABSTRACT 

With the forthcoming International GPS Service (IGS) campaign scheduled for June 21, 
1992 to September 21, 1992, and with the expectation that Analysis Centers will provide products 
within two weeks of the collected data, an examination of strategies for the generation of potential 
products is appropriate. The primary products of the IGS arc GPS ephemerides. Earth rotation, 
baselines and reference frame information. This paper describes an analysis of GIG-9 1 data for the 
purpose of examining strategies for the generation of IGS products. Preliminary results from the 
studies has shown agreement with Lageos Earth rotation results at the 0.8 mas and 1.0 mas levels 
m pole position (x,y), respectively. Other strategies have produced results at the 1.5 mas level. 
Comparison of selected baselines with those obtained in other campaigns has shown agreement at 
the level of 10 ppb, and agreement with baselines determined by VLBI at the several ppb level. 


1. INTRODUCTION 

In January and February, 1991, one of the most ambitious global GPS campaigns to date was 
undertaken, known as GIG-91. This campaign [Melbourne, 1992] included a variety of receivers 
in most areas of the world. For the first time, almost 20 high quality Rogue P-code receivers were 
used at global sites. In some sense, the GIG-91 was a precursor for the 1992 IGS Campaign, which 
will commence on June 21 and end on September 21. An additional campaign, known as EPOCH 
‘92, centered on August 1 , will provide an opportunity for a variety of regional activities. The IGS 
concept is described by Mueller and Beutler[1992]. 

Potential products of the IGS have been extensively discussed and the report of a panel 
charged to identify those products and the timely availability is given by Schutz et al. [1991]. In 
summary, the panel noted that Earth rotation, GPS ephemerides and reference frame/baselines 
would be products with the widest utility. Timeliness of the products was deemed important and 
the expectation that some products could be available within a few days to a few weeks was noted. 

With this background, the primary purpose of this paper was to conduct experiments using 
the GIG-91 data set to evaluate estimation strategies that could be used in the IGS. An additional 
purpose was the comparison of baseline results with those obtained by other techniques and from 
other GPS campaigns as a means of assessing the accuracy. 



2. ANALYSIS SOFTWARE 


All software used in the analysis of the SWP-90 data has been developed at the Center for 
Space Research (CSR) and is known collectively as TEXGAP (TEXas Gps Analysis Programs). 
The analysis process is divided into a preprocessing component and a geodetic component. In the 
preprocessing component, the data were reviewed and corrected for cycle slips, erroneous points 
and general data anomalies. In this process, the time tags of the phase measurements were validated 
and/or corrected using the Lj C/A pseudo-range, or Lj/L^ if the receiver operates with the P-code. 
Finally, explicit double difference ionospherically-corrected measurements were formed for the 
geodetic processing stage. 

The geodetic processing was performed using MSODPl(Muld-SatelIite Orbit Determination 
Program), In the general application of MSODP 1 , the GPS epoch orbit elements and selected force 
model parameters were simultaneously estimated with three-dimensional coordinates of the GPS. 
This software has undergone comparison with programs used for precision orbit determination of 
geodetic satellites, such as Lageos, Starlette and Etalon, all of which are targets for precision sat- 
ellite laser ranging instrumentation [Tapley et al., 1985]. 

3. DATA AND MODELS 

As previously noted, the GIG-9 1 data were used for the study as shown in Table 1 . Although 
the network is dominated by Rogue receivers, TI-4100 receivers at W. Samoa and Easter Island 
were included to improve the southern hemisphere and Pacific coverage. In addition, a Minimac 
2816 at Hobart was included because of the availability of a survey tie to VLBI at the time the in- 
vestigation began; however, the local surveys have recently become available for Yaragadee and 
Tidbinbilla/Canberra. 

The GPS force models followed the current IERS Standards [McCarthy, 1989] and scale fac- 
tors on the ROCK4 [Fliegel, 1992] and y-bias parameters were estimated. The Chao [1973] tropo- 
sphere model was used, and zenith delay parameters were estimated at 2.5 hour intervals from all 
sites. Lageos-derived polar motion and UT1 were used as a priori. For all cases, dual frequency 
double differenced phase measurements were used in the analysis. Pseudo-range measurements 
were used to verify and/or correct the respective receiver clocks. 

The reference frame was based on Lageos satellite laser ranging (SLR) analysis, CSR91L03 
[Eanes et al., 1991] and Very Long Baseline Interferometry (VLBI) analysis GLB718 [Ma et al., 
1991]. The VLBI sites were transformed into the SLR reference frame using transformation pa- 
rameters derived from 18 common sites. The technique has been described by Ray et al.[1991], 

4. ESTIMATION STRATEGIES 

For this study, two primary estimation strategies have been used and a third strategy was par- 
tially examined. The first strategy was based on independent one-day arcs in which orbit parame- 
ters (including y-bias and solar radiation pressure parameters). Earth rotation parameters (x,y) and 
station coordinates were estimated without a priori constraints (i.e., the a priori covariance was es- 
sentially infinite). The second strategy was based on a multi-day estimation of station coordinates. 



but daily determinations of orbit, force model and Earth rotation parameters. For some cases that 
used this strategy, a priori covariance constraints were used. The third strategy was based on a 
multi-day orbital arc, but daily solutions for station coordinates and Earth rotation parameters were 
obtained. 

It is well-known that the models used to describe the dynamics of the GPS satellites are in- 
complete or contain errors (or both). These model deficiencies will lead to a discontinuity at the 
common time point between successive one-day arcs and will produce higher RMS measurement 
residuals on multi-day arcs unless the model deficiency is accommodated by estimated parameters. 
The latter accommodation of errors may not produce improved model parameters as the unmod- 
eled effect may have a signature similar to the other effects, thus allowing the model error to be 
absorbed in other parameters. 

In the multi-day station coordinate strategy, all stations were allowed to adjust, also referred 
to as a “free fiducial” case by Blewitt et al.[1992] and others. This strategy leads to a very ill-con- 
ditioned, or nearly singular, problem when Earth rotation parameters are estimated also and re- 
quires the introduction of some a priori constraints. The constraint commonly used is an a priori 
covariance with coordinate uncertainties chosen to be a specified value, e.g., 100 meters. Other 
ways of avoiding the singularity are to fix the coordinates of some stations or a combination of co- 
ordinates at more than one station. The minimal number of constraints required depends on the pa- 
rameters being estimated. 

5. EARTH ROTATION RESULTS 

Using the first strategy of independent one-day arcs and fixing the coordinates of Hobart, 
Goldstone and Wettzell to the values given in Table 2, the GPS orbit parameters, other station co- 
ordinates and (x,y) Earth rotation were estimated. Although all of the fixed sites were at VLBI lo- 
cations, Hobart is a Minimac receiver and the other two sites use Rogue receivers. In any case, the 
RMS differences of the estimated rotation pole position, compared to Lageos values [Eanes et al., 
1991] produced RMS differences of 1.5 mas in x and 1.4 mas in y after removal of a 5 mas bias. 

In an alternative case, Yaragadee, Goldstone and Wettzell were fixed and all non-Rogue sites 
were eliminated. The RMS differences in pole position were 2.6 mas in x and 3.3 mas in y. Further 
investigation is required to determine whether the cause of the change is associated with the fixed 
coordinates of Yaragadee or with the exclusion of the non-Rogue receivers. 

Usin S strategy in which a multi-day solution was obtained for the stations in a “free fidu- 

cial mode with 100 m a priori on the station coordinate covariance elements, 1 day solutions for 
pole position (x,y) were obtained using 10 mas a priori. The RMS differences of the pole position, 
compared to Lageos, were 0.8 mas in x and 1 .0 mas in y. In all comparisons, the RMS differences 
will change slightly if the GPS results are compared against other Earth rotation series. 

6. BASELINES 

From the three fixed site strategy, the daily repeatability for baseline length on selected base- 
lines is shown in Table 3. The selected cases are all cases in which double differences were directly 



formed for the solution process. It can be noted that the repeatability for all cases involving the TI- 
4100 receivers exhibit worse repeatability than the other cases. It should be noted that for both W. 
Samoa and Easter I. the preprocessing identified some significant systematic features that could be 
related to multi-path problems. 

An additional case was examined in which Trimble receivers at Wellington, New Zealand, 
and Townsville, Australia were included for the puipose of estimating the coordinates of these 
sites. In the case of Wellington, a result obtained during a 1990 campaign afforded an additional 
comparison. Data collected during July 1990 and processed as part of the Southwest Pacific Project 
(SWP) provided a set of coordinates for Wellington. The SWP results [Schutz et al., 1992] used a 
global network that differed from the GIG-9 1 in two primary ways. First, the SWP global network 
was dominated by “codeless” receivers and, second. Selective Availability (S A) was activated. As 
noted previously, the GIG-91 global network was dominated by Rogue P-code receivers and SA 
was not implemented. The comparison of the Hobart to Wellington baseline is given in Table 4 and 
the coordinates of the Australia/New Zealand sites derived from GIG-91 are given in Table 5. 

7. CONCLUSIONS 

Based on the preliminary results given in this paper, it has been shown that Earth rotation 
components (x,y) were obtained that agree with other determinations at the 1 mas level (RMS). 
Two strategies were examined: a multi-day case and cases using independent one-day arcs. Base- 
lines from the three fixed site case show repeatability at the several ppb level, except for cases us- 
ing the TI receivers in the Pacific which are at the level of 10-20 ppb. Further examination of the 
influence of mixing Rogue and TI data will be conducted and other estimation strategies are under 
examination. 

Results for Wellington that were obtained from two campaigns show agreement at the 10 ppb 
level. In one case, the campaign was dominated by global codeless receivers, whereas the GIG-91 
was dominated by Rogue receivers. An additional difference was the fact that SA was activated 
during the earlier campaign, but not during GIG-9 1 . 
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TABLE 1. GIG -91 SITES 


Rogue Receivers: 

Yaragadee, Australia 
Canberra, Australia 
Santiago, Chile 
Hartesbeestoek, S. Africa 
Kokee Park, Hawaii, 

Usuda, Japan 
Goldstone, CA 
Victoria, BC 
Fairbanks, Alaska 

Minimac 28 16 AT Receiver 

Hobart, Tasmania, Australia 

TI-4100 Receivers: 

W. Samoa 
Easter Island 


Wettzell, Germany 
Madrid, Spain 
Matera, Italy 
Kootwijk, Netherlands 
Ny Alesund, Norway 
Tromso, Norway 
Algonquin, Ontario 
Yellowknife, NWT 


TABLE 2. COORDINATES USED FOR FIXED SITES (m) 


Site 


y 


Goldstone (Rogue -2353613.9840 

Hobart (Minimac LI) -3950 1 84.0724 
Wettzell (Rogue) 4075579.3868 


-4641385.4730 

2522364.5271 

931807.2475 


3676976.4990 

-4311588.6675 

4801570.9395 


Note: The Rogue coordinates refer to the top of the antenna 


TABLE 3. SELECTED BASELINES 


Sites 

Baseline Length 
(km) 

Daily RMS Scatter 
(ppb of baseline length) 

Algonquin - Wettzell 

6154 

5.6 

Goldstone - Yellowknife 

2986 

8.8 

Goldstone - Fairbanks 

3807 

6.3 

Goldstone - Algonquin 

3402 

12.4 

Yaragadee - Tidbinbilla 

3197 

8.9 

Tidbinbilla - W. Samoa 

4474 

17.7 

Kokee Park - W. Samoa 

4124 

22.2 

Kokee Park - Usuda 

5894 

8.4 



Case 


TABLE 4. COMPARISON OF HOBART-WELLINGTON BASELINE (m) 


x 


y 


L 


SWP-90 

GIG-91 


-830464.717 

-830464.699 


-2085857.511 

-2085857.541 


126148.175 

126148.187 


2248641.049 

2248641.070 


Difference (xyz, L) 0.018 -0.030 0.012 

Difference (NEU) 0.005 0.028 0.023 


0.021 


(NEU: North, East, Up) 


TABLE 5. COORDINATES OF SELECTED AUSTRALIA/ NEW ZEALAND SITES (m) 
Data: GIG-91 


Site 

X 

y 

z 


Hobart 

-3950184.072 

2522364.527 

-4311588.668 

Minimac L 

Wellington 

RMS 

-4780648.771 

0.051 

436506.986 

0.017 

-4185440.481 

0.032 

Trimble L 3 

Tidbinbilla 

RMS 

-4460987.995 

0.025 

2682362.260 

0.013 

-3674626.550 

0.033 

Rogue L 3 

Townsville* 

RMS 

-5041024.956 

0.020 

3296980.304 

0.039 

-2090553.463 

0.009 

Trimble L 3 


RMS refers to the daily scatter in the solutions 
* denotes that some solutions were edited 
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GPS Reference Frames and Orbit Accuracy 

B E Schutz . PAM Abusali, M Watkins, H Rim, D Kuang (Center 
for Space Research, The University of Texas at Austin, Austin, TX 
78712; ph. 512-471-5573; fax. 512-471-3570) 

Analysis of GIG-91 data has been used to investigate nongravitational 

forces acting on the GPS satellites and their influence on the 
determination of geodetic parameters. It has been previously demonstrated 
that significant differences exist between the eclipsing and non-eclipsing 
satellites, which have prompted examination of eclipsing phenomena 
and associated forces. The analysis has been aided by the extended 
global tracking network available during the GIG campaign. The 
influence of the orbit accuracy on the reference frames is also examined, 
and comparisons between common sites with other campaigns will be 
made, using baseline lengths ranging from a few hundred kilometers 
to several thousand kilometers, with attention given to the determination 
of the vertical component. 




